/ 9 ? ~T /!/ ¥ f / 

IMPROVEMENTS AND NEW FEATURES IN THE PDF MODULE 
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Cleveland, Ohio 


Overview 


• Modeling: What models are used in this package and what 
are their advantages and disadvantages. 

• Numerics: Describe how the PDF model is implemented 
and what are the features of the program. 

• Future Developments: What can be expected in the future 
from the NASA Lewis PDF code. 
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PDF Modeled Equations. 


• Exact scalar PDF transport equation is: 
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• Terms on the LHS exact - need to model the four terms 
on RHS, corresponding to turbulent convection, molecular 
mixing and the pressure term. 


Turbulent Convection 


• This term is modeled as a simple gradient diffusion process. 

(pUi\±,v)P ~ ^ 

• Dt is the turbulent diffusion coefficient, equal to the eddy 
viscosity. (Assume unity Schmidt) 

• Disadvantage: Counter-gradient diffusion known to occur 
in some pre-mixed flames. 


270 


Molecular Mixing. 


• Molecular mixing can be viewed as process which changes 
the shape of the scalar PDF without affecting the mean. 

• Molecular mixing is modeled by two models: A coalescence/ 
dispersion model (Hsu and Chen) and a relax-to-mean model 
(Dopazo). 

• Advantages of both models is that they are simple and 
readily adaptable to any number of scalars. 

• Disadvantages are the relative lack of physics in the models. 


Pressure Term. 


• Pressure term model is based on second order closure mod- 
els for compressible flows (eg. Sarkar). 
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• Advantages are that model is tried and tested in finite 
volume codes. Disadvantage is that only the mean pressure 
can be used for model. Idealy we would like a stochastic 
process for two state variables. 
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Numerics 


• Solution of scalar PDF transport equation achieved by a 
particle based Monte Carlo scheme. 

• PDF represented by an ensemble of particles, each with a 
composition and enthalpy. 

• PDF evolves by the motion of these particles in physical, 
scalar and enthalpy space, by exact and modeled processes, 
eg. Convection, reaction, mixing. 

• Statistics (eg. means) obtained by averaging over ensemble 
of particles. 


Numerical Details - Monte Carlo Scheme 

• Module based on cell-centered quantities. 

• PDF method is a nodal one. ie. All particles reside at the 
center of the cell, and can move only to neighboring cells. 

• Number of particles at each node is the same and the 
number remains fixed. 

• Statistics obtained by averaging over particles at each node, 
and also by time-averaging. 
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Numerics - Convection 


• Upwind scheme used for mean convection, and central dif- 
ference for turbulent convection. 

• Evolution achieved by moving particles in from adjacent 
nodes. Particles are selected at random. 

• Fractions of particles are treated by random convection. 


if 6.3 particles then 


6.0 70% of the time 

7.0 30% of the time 


(4) 


Numerics - Reaction 

• Although reaction source term treated exactly, several dif- 
ferent numerical schemes are needed. Timing figures are for 
% of time spent in the PDF part of the code on SPARC II 
workstation. 

• No reaction: For scalar mixing calculations. 

(Timing: 17.2 %) 

• Equilibrium reaction: Assume reaction proceeds at infinite 
speed. Table of equilibrium composition as a function of 
mixture fraction obtained from seperate CHEMKIN routine. 
(Timing: 24.4 %) 
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Numerics - Reaction Cont. 


• One-step global reaction schemes. Westbrook and Dryer 
global reactions integrated for each time-step. 

(Timing: 51.1 %) 

• Tabulated reaction incriments. Users create their own table 
of composition incriments as a function of scalars using the 
adaptive tabulation scheme provided, plus the users favourite 
reduced mechanism. 

(Timing: 58.9 %) 

• Chemkin full mechanism integration. Very slow and not 
recomended except for parallel applications. 

(Timing: 97.8 %) 


Numerics - Averaging 

• To reduce statistical error in the evaluation of the mean 
scalar quantities (without increasing the number of particles 
per node), time averaging is employed. 

• A weighted time average is used to give more weight to 
recient values and less to those in the far past. 

(( <p))t = 3— r((4>)t + ( 5 ) 

Wt T 1 

wt = ct(wt - 1 + 1 ) ( 6 ) 
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Numerics: Misc. 


• A portable random number generator is now included in the 
module, set up for 32 bit machines. 

• A time step check is now performed to ensure boundedness 
of the PDF solution, ie. no negative numbers of particles. 

• Rplus/PDF release ported to workstation enviroment. K- 
epsilon now standard turbulence model. 


Future Work. 

• Release of 3D version with new improvements. 

• Implimentation of parallel processing for distributed cluster 
environment. (PVM based) 

• Include model for another state variable to close PDF mod- 
eling. 
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